Neutrinoless double beta decay at the LHC 
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We compare the discovery potential of the LHC for lepton number violating (LNV) signals with the 
sensitivity of current and future double beta decay experiments, assuming 0vf3p decay is dominated 
by heavy particle exchange. We consider charged scalar, leptoquark and diquark mechanisms of 
OvPfl decay, covering the Ou/3/3 decay operators with both, the smallest and largest, possible rates. 
We demonstrate, if Ou/3/3 decay were found with a half-life below 10 26 - 10 27 ys a positive signal 
should show up at the LHC, except for some particular cases of the leptoquark mechanism, and vice 
versa, if the LHC does not find any hints for LNV, a "short-range" explanation for a finite Ou/3/3 decay 
half- life will be ruled out in most cases. We argue, if a positive LNV signal were found at the LHC, it 
is possible to identify the dominant contribution to 0vf3/3. Two different kinds of observables which 
could provide such "model discriminating" power are discussed: Different invariant mass peaks and 
the charge asymmetry. 



1. Introduction: Neutrinoless double beta (Ozz/3,3) de- 
cay can be induced by the exchange of light particles 
such as Majorana neutrinos (long-range mechanism) or 
by heavy particles (short-range mechanism). The clas- 
sical example of the latter is right-handed W-boson and 
heavy neutrino exchange, which occurs in left-right sym- 
metric extensions of the standard model [H]. A number 
of other examples of such "short-range" contributions to 
Qvj3p decay have been discussed in the literature (for a 
review see |2j ) , but only very recently the general decom- 
position for the Qvj3j3 decay operator, O^J 1 ^ , in terms of 
rcnormalizablc vertices connected by heavy virtual parti- 
cles, has been given in Q. Interestingly, several new dia- 
grams for 0vf3j3 decay were identified in this work, many 
of which require quite exotic particles such as diquarks 
or coloured fermions with fractional charges. Common 
to all of the short-range contributions is that masses of 
order M ~ (1 — 5) TeV are needed for a 0v[3f3 decay 
half-life in the range of T°J^ ~ (10 25 - 10 27 ) ys in either 
76 Ge or i36 Xe _ 



The very same diagrams which lead to Ov/3/3 decay can 
lead to LHC signals with like-sign dilcptons plus (at least) 
two jets in the final state HQ. The ATLAS § and the 
CMS collaborations have published searches for such 
events. (Disregarding flavor, we will refer to these events 
as "eej j-like" .) Assuming on-shell production of a right- 
handed W-boson, Wr, decaying to a heavy neutrino, N, 
plus lepton they then interpret their null result as 
lower limits on m WR of m WR > (2.5 - 2.9) TeV Q, de- 
pending on mjv and assuming the coupling of Wr 
to fermions is equal to the standard model coupling gr,. 



Note that provides separate limits for opposite-sign 
(OS) and same-sign (SS) dileptons, while CMS provides 
only combined limits. While, of course, the search for 
OS and SS is experimentally very similar, obviously only 
SS type events are a proof for the discovery of lepton 
number violation (LNV). 

In this paper we discuss how searches for LNV signals 
at the LHC compare with the sensitivities of current and 
future 0vfif3 decay experiments. We do not limit our- 
selves to a particular model of LNV, but use a general 
decomposition of 0^=9 ■ ^ non-observation of eejj-like 
events in the future LHC run at = 14 TeV would 
rule out all but the pure leptoquark decompositions as 
the dominant contribution to 0^=9^ ■ Perhaps even more 
interesting is that, if a positive LNV signal were found at 
the LHC, it is possible to identify the dominant contribu- 
tion to Off3(3. We discus two types of observables which 
could provide such "model discriminating" power: The 
construction of different invariant mass combinations and 
the ratio of positron-like to electron-like events. We will 
call the latter the "charge asymmetry" . 

2. Short-range double beta decay: Double beta decay 
can be induced by exchange of scalars or vectors with 
or without an accompanying virtual fermion. Since re- 
sults for scalars and vectors are qualitatively very similar, 
we will concentrate on the scalar case only. Also, there 
are two basic topologies of the diagrams contributing to 
Ovpp decay. Here, we limit ourselves to the discussion 
of topology-I (scalar-fermion-scalar exchange) shown in 
fig. ([!])■ Results for topology-II (triple-scalar exchange 
coupled to the trilinear vertex) will be presented else- 
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where [8|. There are in total 18 different possible decom- 
positions ( "diagrams" ) in topology-I that can contribute 
to the 0vpf3 decay operator, O^Jg . Depending on the 
color and SU(2)l transformation properties of the in- 
ternal fermion and bosons, each decomposition can be 
realized in several ways, see below. 12 decompositions 
involve fractionally charged (and colored) fermions, and 
were discussed for the first time in The remaining 6 
possibilities have all been discussed before in the litera- 
ture in the context of concrete models, most notably the 
mass mechanism (for a recent review see Q) or R-parity 
violating supersymmetry If], 11 1. 
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FIG. 1: Example diagrams for short-range double beta decay, 
see text. 

In fig. (dJ we show some example diagrams for 0isf3f3 
decay. These examples were chosen such that all six dif- 
ferent bosons (S+i, S+2, £2/3 > ^1/31 ^4/3 an< ^ ^2/3) anc ^ 
all four types of fermions (-00, ^1/3, 1P4/3 and ^5/3), which 
can contribute to CV/3/3 decay, appear at least once each 
in the figure. Note that replacing S+i in diagram (a) on 
both sides by a standard model W-boson and interpreting 
ipo as neutrino would lead to the classical mass mecha- 
nism and diagrams of the type shown in (b) appear, for 
example, in R-parity violating SUSY, while diagrams (c) 
and (d) are examples of the new type of diagrams. 

Integrating out the heavy particles in all cases leads to 
an amplitude 
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We denote the couplings as g\ — (74, arbitrarily read from 
left-to-right in the diagrams. All of these are in principle 
different free parameters in general, mg^ and rn^, k 
are the masses of the bosons and the fermion, respec- 
tively, not necessarily mg i = msj ■ 

While 0^=9^ depends on the particle physics param- 
eters always in the form of eq.(JT]), the nuclear matrix 
elements (and thus the half-life) depend on the chirality 
of the "outer" fermions as well as the colour transforma- 
tion properties of the "inner" particles in the diagram. 



The recent paper Q provides a complete list of oper- 
ators (for short-range Qvf3j3 decay with scalars), which, 
once numerical values for the nuclear matrix elements are 
taken from, for example [liq . can be directly converted 
into calculated half-lifes. Current limits from 76 Ge 
and 136 Xe 0, 



15|, both of the order of T"^ p > 10 
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TeV, 



ys, correspond to roughly M e ff > (1.2 — i.2)g e ^ 
depending on decomposition. 

3. Sensitivity: LHC versus 0v(3f3 decay: The number 
of eejj-like events at the LHC in general depends on a 
different combination of couplings and masses than the 
0^/3/3 decay amplitude. Here, for brevity, we will discuss 
only the "symmetric decompositions", see below. We 
stress, however, that we have checked that, despite being 
a special case, the results presented below still cover both, 
the most optimistic and the most pessimistic, scenario 
for the LHC. A systematic case- by-case study covering 
all decompositions and cases in detail will be presented 
elsewhere 

In "symmetric" decompositions, such as diagram (a) 
of fig. (p}, couplings are pairwise equal, i.e. g\ = 174 and 
g2 = 93- Assuming that such an equality holds at least 
approximately, in the small width approximation we can 
write the number of events at the LHC as 



#(eejj)/£ = a(pp Si) x Br(S; -)■ eejj) (2) 
= F Si 



m^ k 



liffBv^iSi^eejj). 



Here, we have rewritten Fg i (ms i ) = cr(pp — > Si)/g\ 
and Br e// (S' i -> eejj) = Br(5< ->• eejj)/g%. While 
Br(5j — ► eejj) is fixed for a given decomposition, once 
all couplings and masses are fixed, here we simply 
treat Br 6 " (Si —> eejj) as a free number smaller than 
one. Note, however, that in the limit where all cou- 
plings are equal (and m^ k = 0) Br(5+i — > e + e + jj)= 
Br(S_|_i — > e + e~ jj) ~ 1/8 for S+±. In the LQ case, the 
LQ is produced in association with a lepton, i.e. we cal- 
culate a(pp — > Si + e) x Br(Si —> ejjj), see next section. 
We have used CalcHEP and MadGraph 5 to cal- 
culate the production cross sections for S+i, S LQ , and 
S D ® at the LHC. We have compared our results with 
the literature 
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- 20j | and found quite good agreement in 
all cases. 

In fig. ([2]) we then plot the compared sensitivities of 
0^/3/3 decay and the LHC for five different cases, using 
eq. @. The more complicated "asymmetric" case, with 
more than one scalar mass and all couplings different, will 
be presented elsewhere [|[ , but leads to the same conclu- 
sions. For the LHC we show the expected sensitivity lim- 
its for a x Br(eejj) = 1 fb (dashed lines) and 0.1 fb (solid 
lines) in the plane <? e // versus M e ff, for two different val- 
ues of = 200 GeV (top) and = 1 TeV (bottom) . 
The different color codes correspond to the five different 
scalar bosons, than can be singly produced at the LHC, 
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FIG. 2: Sensitivity limits for the LHC for production of five 
different scalar bosons compared to current and future double 
beta decay experiments, see text. 

namely S+i, Sy®, S^®, and Sy® ■ Note that S+i is 
very similar to the well-known case of Wr production in 
LR symmetry and S^H and Sf® correspond to the most 
pessimistic and most optimistic scenario for the LHC. 

In addition fig. <j2j) shows four different cases for 
current and future limits from Ov/3/3 decay. The dark 
gray area is the currently excluded part of parameter 
space from non-observation of 136 Xe decay with > 
1.6 x 10 25 ys [3| assuming 0^/?/3 decay is caused by the 
operator with the smallest rate (j^Oi in the complete 
list), and thus corresponds to the most pessimistic case 
for the sensitivity of 0v/3/3 decay. The three blue areas 
are (from left to right): Smallest rate, but for a limit 
of r°" 3 w > 10 26 ys, largest rate with ifff? > 10 26 ys 

and, finally rate with > 10 27 Y s - The lightest area 

to the right therefore corresponds to the most optimistic 
reach for 0vf3/3 decay in the foreseeable future. 

As can be seen, with the exception of the LQ cases, 
the LHC at = 14 TeV will be more sensitive than 
(V/3/3 decay experiments as probe for LNV. For the LQ 
case, the LHC is more sensitive than Gv/3/3 decay in the 
pessimistic case for Gv/3/3 (operators 0\ and O5 in the 
notation of [l2[) but not for the one to which Gv/3/3 de- 



cay is most sensitive, particularly O^- For the remaining 
operators O3 and O4 CV/3/3 decay and LHC sensitivities 
are very similar. Note that decompositions with one LQ 
and one S + i or one S* 3 ® will be similarly constrained by 
LHC searches as the "symmetric" 5+i or Sf® case. 

4- Distinguishing LNV models: As long as LHC ex- 
periments provide only upper limits on LNV signals, 
conversion to model dependent exclusion plots is, apart 
from different numerical factors, quite similar in all cases. 
However, if a positive signal is observed in the future, the 
LHC has the potential to differentiate between different 
LNV models. The two probably most promising observ- 
ables to do so arc: (i) different invariant mass peaks and 
(ii) the charge asymmetry. 

4 .a Invariant mass peaks: There are three different 
sub-cases to discuss: (i) Decompositions with either <S-|_i 
or Sf®; (ii) decompositions with Sf®; and (hi) decom- 
positions with fractionally charged fermions. 

Consider first the simplest case, diagram (a) of fig. 
(Q]). This case is equal to the well-known LR-symmetric 
case: On-shcll production of S + i decaying to -00 + e , 
with tp ~ > ejj via off-shell (S+i)*, leads to two invari- 
ant mass peaks: to 2 , = m% and m 2 •_■ ~ m 2 ,, . On 
the other hand, m 2 e shows a broad distribution. Com- 
pare this to, for example, fig- (jT|) diagram (c). Here, still 
m lejj = TO s +1 ' but m\ en = to 2 o and since ms +2 < m x p i 
is possible, a third peak can show up at m 2 e = rng +2 . In 

general, all decompositions with S+i or S["^ will have 
a mass peak at ni/^ e jj and different combinations of sub- 
systems will also have peaks, depending on the exact de- 
composition and mass ordering. 

The case (ii), decompositions with LQs is different, 
since the LQ is produced in association with a lepton, 
i.e. there are (at least) three jets in the final state and 
m-eejj does not show a peak at any mass (neither does 
TO Ljjj)' wm m l 2 jjj = m2 S LQ- Again, different subsys- 
tems of e2jj j show peaks, depending on mass ordering 
and exact decomposition, but in the LQ case there never 
is a peak in the m 2 e distribution. Finally, as mentioned 
above, for case (hi), 12 decompositions have fraction- 
ally charged fermions, all of which have to be coloured. 
These fermions can be pair-produced 0) and the final 
state of their decay is eejjjj. A mass peak is found in 
m eejjjj = ^ m %i anc ^' m addition, subsystems such as 
m 2 ^- or m\ e j or m?^ must show peaks at . . 

Note, that the three sub-cases are not mutually exclu- 
sive: A decomposition with a diquark and a leptoquark 
with different masses, for example, could show multiple 
peaks in the spectrum. Also note that all diquark de- 
compositions (and many others) do come with coloured 
fermions, thus eejj-like events can be accompanied by 
eejjjj-like events with peaks at in general different in- 
variant masses in many cases. 

4-b Charge asymmetry: Define the charge asymmetry 
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FIG. 3: Charge asymmetry, i.e. ratio of positron-like to 
electron-like events as a function of the boson mass for differ- 
ent kinds of scalars. 



TeV, but realistically (compare fig. ([2])) for LQs the LHC 
will not have sufficient sensitivity to find any events for 
LQ masses larger than roughly 2 TeV. 

5. Conclusions: We have compared the sensitivity of 
the LHC experiments for LNV signals with current and 
future sensitivities of double beta decay. The LHC has 
large discovery potential for general LNV violating mod- 
els. In case a positive signal is found in the future at the 
LHC, we have proposed two different observables, which 
allow to distinguish different LNV models: (i) construc- 
tion of different invariants mass peaks and (ii) the ratio 
of positron to electron same-sign dilepton events. 
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as the ratio of the number of events 

xca = #(e + e+)/#(e-e-). 



(3) 
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Consider first the simplest case of diagram (a) in fig.fJTJ. 
With only one boson and two pair-wise equal couplings 
(.9i = .94 an d .92 = .93) xca is simply determined by 
the ratio of production cross section for o~(pp — > S+i) 
divided by cr(pp — > S-i). Since the number of u quarks 
is different from the number of d quarks in the proton, 
the two cross sections are not the same and one expects 
to see more events with e + e + than e~e~. More generally, 
xca is a function of all masses and couplings, but even in 
the limit where couplings are equal, different Ov/3/3 decay 
decompositions have different values for xca- 

Fig. ([3]) shows the theoretically expected value of xca 
for six different cases in the equal coupling limit. These 
cases shown do not (always) directly correspond to a 
particular decomposition, but represent extreme limits, 
where the scalar identified by the colour code gives the 
dominant contribution to xca- The case S®® + S®® as_ 
sumcs that both diquarks have the same mass and cou- 
plings and thus is just the average of the two individual 
cases. If both 5^5 and contribute equally to the 
Ov(3(3 decay rate (case not shown), the two asymmetries 
average to a number close to, but not equal to the case of 
5+1. Note that the case for S+i is equal to the expecta- 
tions in the left-right symmetric model. In this cases, the 
charge asymmetry varies only weakly from xca — (2 — 4) 
in the mass range shown. However, as the figure demon- 
strates, in other cases much larger (and much smaller) 
xca are possible. Experimentalists should therefore pro- 
vide data sets separately for positively and negatively 
charged leptons, in order to cover all LNV models. Fi- 
nally, note that for the LQ case we show curves up to 6 
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